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We present a characteristic algorithm for computing the perturbations of a 
Schwarzschild spacetime by means of solving the Teukolsky equations. Our meth- 
ods and results are expected to have direct bearing on the study ofrbinary black 
holes presently underway using a fully nonlinear characteristic code El. 

Recent years have shown a vital synergism between numerical and perturbative 
approaches in the context of Cauchy evolution. In this work we try to make these 
benefits also accessible in the context of null-surface-based evolution. We treat the 
"close-limit" approximation to the post-merger phase of a binary black hole space- 
time in terms of the characteristic, initial value problem. We outline a framework 
to treat the Teukolsky equation in a null evolution setting, and our numerical 
strategy to solve the perturbative problem in the time domain.pA detailed discus- 
sion and results for the outgoing radiation problem are given in a, the details of the 
ingoing radiation problem will be treated in a forthcoming paper □. 

We apply our method to evolve "close-limit" horizon data, which we obtain from 
linearizing the data discussed in □. There we use Bondi-Saschs coordinates, which 
simplify to Israel coordinates in the spherical case, which are based on an afhne 
parameter u along the white hole horizon, normalized to u ~ at the bifurcation 
sphere, and an afhne parameter A along the outgoing null cones, fixed to A = 
at the white hole horizon and (d a u)(d a \) = —1). These coordinates cover the 
entire Kruskal manifold r > with remarkably simple analytic behavior, as first 
discovered by Israel 13. Using these convenient coordinates as a starting point, it 
is straightforward to write the Teukolsky equations explicitly in any background 
coordinate system. 

In order to treat the radiation near X + it is advantageous to consider a "boosted 
tetrad" (l a , h a , m a , fh a ) with l a = — V a u, n a satisfying l a n a = —1 and a dyad 
171(1 = v^t'-m) + snh?(^) a - W e accordingly define boosted Weyl scalars ipo — 
C a bcdl a 'm b l c 'm d and yj^ = C a bcdn a m h h c m d . This boosted tetrad is adapted to the 
afhne time on u at I + rather than the afhne time u at the horizon. 
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This tetrad has proven useful for the evolution of outgoing radiation described by 
ip4, for ingoing radiation described by ipg, a more sophisticated choice is necessary 
to achieve optimal results, as will be discussed in □. 

In order to accurately track the evolution of ipi well into the power law tail phase 
(which typically requires more than a 1000 M) of evolution in Bondi time u, we 
factor out the asymptotic falloff of ip4 for large r, we use fixed time steps in Bondi 
time and a compactihed version p of the tortoise coordinate r* = r + 2M \n(r /2M — 
1), p = potanp, where po is an adjustable parameter. The coordiante p allows for 
a fixed grid to accurately resolve the Schwarzschild potential at all times, which 
proved vital in the tail phase. For a typical choice of initial data the power law 
tail only sets in after the quasinormal oscillations have decayed by more than 10 
orders of magnitude. In order for the final tail not to be lost in machine error it is 
necessary to evolve the quasinormal phase in quadruple precision. 

We have obtained preliminary results to construct close approximation black 
hole null data without ingoing radiation, by using superposition of solutions for the 
linear Teukolsky equation, and by utilizing the time symmetry of the background 
Schwarzschild spacetime which implies a correspondence map of a retarded solution 
(no incoming radiation) for into an advanced solution (no outgoing radiation) 
for ^o- Details will appear in □. 



Acknowledgments 

M.C. was partially supported by a Marie-Curie Fellowship (HPMF-CT-1999-00334) 
This work has been partially supported by NSF PHY 9988663 and NSF PHY 
9800731 to the University of Pittsburgh. Computer time for this project was pro- 
vided by the Pittsburgh Supercomputing Center and by NPACI. 



References 



1. R. Gomez et. al., Phys. Rev. Lett. 80, 3915-3918 (1998). 

2. S. A. Teukolsky, Astrophys. J. 185, 635-647 (1973). 

3. Manuela Campanelli, Roberto Gomez, Sascha Husa, Jeffrey Winicour and 
Yosef Zlochower, |gr-qc/0012107 . 

4. Manuela Campanelli, Roberto Gomez, Sascha Husa, Jeffrey Winicour and 
Yosef Zlochower, to be published. 

5. Roberto Gomez, Sascha Husa and Jeffrey Winicour, gr-qc/0009092j . 

6. Werner Israel, Phys. Rev. 143, 1016 (1966). 



TeukTalk: submitted to World Scientific on February 7, 2008 2 



